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contributed approximately 85.81 % of the overall exergy destruction. The 
system's exergo-thermal index (ETI) stood at 1.713, with the GT only having an 
ETI of 2.106. Similarly, the exergetic sustainability index (ESI) of the 
multigeneration plant was not greater than 2.04. The exergoeconomic analysis 
shows a low average energy cost from the GT, estimated at 0.836 $/GJ, 
compared to the Kalina subsystem, which stood at 6.53 $/GJ. The 
thermodynamic and cost evaluation of the system demonstrates substantial 
benefits from the plant, which kept the hydrogen production rate at 0.1524 
kg/hr. 


1. Introduction 


The role of energy availability in economic development 
cannot be over-emphasized. Energy remains the fulcrum for 
sustainable development [1]. For the past decades, the 
fundamental processes for energy generation, especially 
electricity, have been from burning fossil fuels in gas turbine 
plants, steam plants, and a combination of the two in 
cogeneration [2, 3]. Several plant configurations have been 
developed and implemented to increase these thermal plants’ 
energy efficiency. Research is increasing to produce different 
products from these systems, which has increased in scope in 


recent times [4, 5]. However, the configuration of the lower 
bottoming cycles and the choice of operating parameters play 
vital roles in developing sustainable products with minimal 
exergy destruction and environmental impact. Among the 
many options for lower bottoming cycles to increase energy 
conversion system products are the organic Rankine cycle 
(ORC) [6-8], Kalina cycle [9, 10], Goswam cycle [11] steam 
turbine cycle [12], and for water electrolyzers [13]. 
Multigeneration plants have been presented in the literature. 
For example, Dev and Attri [12] presented a Multi-plant cycle 
based on ORC, which produces six products: electricity, 
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heating, hot water, hydrogen, cooling, and dry air. The results 
indicate an energetic coefficient of performance (COP) of 
approximately 6096, while that for the exergetic COP was 
estimated at 1096. Similarly, the study in reference [13] 
presented a primary Brayton cycle with several strings of 
bottoming cycles that can provide electrical power, domestic 
water heating, and cooling through refrigeration processes 
embedded in ORC and vapor absorption cycles. The study 
obtained energy and exergy efficiencies of 44.22 and 61.50 96, 
respectively, with an exergo thermal index of 0.675. 
Additionally, other lower cycles are reported to produce 
hydrogen, while some are used for the desalination of water 
[14, 15]. The study in [14] proposed a multi-cycle power plant 
for cooling, heating, and desalting water production. The 
result shows that the water production rate was 0.364 kg/s 
with a cooling rate and power output estimated at 1 MW and 
30.5 MW, respectively. Similarly, Abam et al. [15] designed 
and developed a novel incorporated cycle based on the solid- 
oxide fuel cell-GT for concurrent fresh water, electricity, and 
hydrogen production. The exergetic efficiency was estimated 
at 54.2 96 at refresh water rate production of 90.1 m3/h at a 
unit cost of 32.9 $/GJ. Furthermore, Anvari et al. [16] analyzed 
an innovative solar-based multigeneration plant. The system 
was developed for power generation, desalination, and 
hydrogen production. The thermodynamic assessment 
indicates that the multigeneration system had 23.2 96 as 
energy efficiency and 6.296 for exergy efficiency. Ref [17] 
performed a thermodynamic and economic evaluation of a 
geothermal-based integrated plant which comprises ORC, 
domestic water heater, electrolyzer, and absorption 
refrigeration system to produce hydrogen and electricity, 
heating and cooling. The total efficiencies were estimated at 
34.9896 energy and 49.1796 exergy, whereas the hydrogen 
production rate and unit cost of production stood at 0.052 g/s 
and 5.967 $/kg, respectively. Furthermore, Khalid et al. [18] 
presented a techno-economic valuation of a solar, 
geothermal-driven multigeneration system to produce; 
electricity cooling, hot water, heating, and hydrogen. The 
hydrogen production rate was estimated at 2.7 kg/h while $ 
476,000 and $ 0.089/kWh were the net presents and the 
levelized cost of electricity. Ref [19] proposed a geothermal 
multigeneration plant comprising an electrolyzer, absorption 
cycle, Kalina, and flash cycles for hydrogen and ice 
production. The results show a maximum exergy efficiency of 
26.25%. Multigeneration systems require adding bottoming 
cycles to gas-fired, geothermal, or solar-based topping cycles 
to generate multiple products. Thus, several choices on the 
mode of a bottoming cycle depend on several factors. These 
may include the exiting temperature of the topping cycle flue 
gas, the required efficiency, the quantity of desired products, 
and the economic and environmental impact of the bottoming 
cycle. The Kalina cycle can provide an efficiency gain of 10-50 
96 compared to the conventional thermodynamic cycles with 
equivalent output However, system modification or 
adjustment in the thermodynamic pathways may be 
necessary to maximize the waste heat used in the bottoming 
cycles. The current study introduced a dual heat input to 
improve the thermal heat requirement of the PEM 
electrolyzer from the same energy source. The latter will 
increase the hydrogen production yield since hydrogen 
production yield is related to the quantity of electricity and 
heat in the PEM electrolyzer. An external thermal heat source 
for the electrolyzer may increase system costs and 
maintenance procedures. In the current study, an external 
heat source for the electrolyzer was avoided. The system was 
designed to utilize the dual heat input from the condensate of 
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the two Kalina condensers. The latter innovation improved 
the hydrogen production yield and reduced the 
multigeneration plant's economic cost. In light of this 
development, this research proposes a modified Kalina-based 
multigeneration cycle incorporating biomass gasification for 
reheating in the combined topping cycle. The study thus 
determines the thermo-environmental and thermo-economic 
performance of the modified Kalina-based multigeneration 
cycle. The best operating parameters corresponding to low 
operating costs and environmental impact were equally 
evaluated. 


2. System description 

The operation starts at state (1) in Figure 1, where the air 
is drawn to the low-pressure compressor (LPC), and the 
temperature and pressure are raised to state (2). An 
intercooler reduces the mechanical work required for 
compression from the state (2) to (3). The cooled air is 
recompressed by a high-pressure compressor (HPC) at state 
(3). The compressed air is heated partly by an expanded gas 
from the low-pressure turbine (LPT). Fuel is added to the 
combustion chamber (CC) at constant pressure, raising the 
energy level of the air stream to the high-pressure turbine 
(HPT), which expands to a pressure sufficient to drive both 
the HPC and LPC. The expanded gas is reheated with biomass 
gas for expansion in the LPT to produce shaft work which 
drives an alternator to produce electricity. After partly raising 
the air temperature in a heat exchanger, the expanded gas is 
made to pass through a vapor generator which powers the 
Kalina cycle embedded with a vapor absorption system (VAS). 
Ammonium water and lithium bromide drive the Kalina and 
VAS system, respectively. In the Kalina cycle, after the 
ammonia water solution receives heat, the energy level of the 
mixture increases at the state (17). The separator is provided 
to separate the ammonium water solution into rich and weak 
solutions. The rich solution expands in the turbine producing 
electrical energy. The expanded vapor, still rich, is further 
separated in another separating vessel. The rich part of the 
mixture is condensed and throttled, and evaporated, 
providing cooling in the first part of this cycle. Meanwhile, the 
weak solution already expanded in the turbine and separated 
is throttled to a pressure equal to the weak solution at state 
(24) and mixes with the hot vapor leaving at state (23) after 
being equally throttled. These two streams heat the desorber 
for powering the VAS. The exiting stream at the desorber is 
passed through a heat exchanger before being added up with 
the resulting stream used for refrigeration at state (31). The 
summation produces a stream at state (32), which passes 
through a heat exchanger and condenser and is finally 
pumped to the vapor generator to commence the next cycle. 
The heat obtained from the dual condensation of the 
expanded refrigerants in the Kalina cycle is used to heat the 
water directed for electrolysis in the PEM electrolyzer. Part of 
the electricity generated from the Kalina turbine is utilized for 
the electrolyser's operation to produce hydrogen. 


3. Methodology and system modelling 
3.1 Thermodynamic modelling 

The general energy flow balances in a thermodynamic 
system under steady for the k^ component are presented in 
Eq (1) [20, 21]. 
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Figure 1. Schematic of the multigeneration plant 
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The general exergy balance for a control volume in a steady 
state, neglecting potential, kinetic, and electrical energy, is 
defined as: 


Exa = Ye (1 = 2) Qk m Wey F Xi(ni£xi) m Ye(neExe) (2) 


"M . TYA ; 
Where E, is the exergy destruction rate, (1 = A) Qp is the 
k 


exergy flow rate accompanying heat transfer, W., is the rate 
of work done within the control volume, n;Ex; and n,Ex, is 
the exergy flow rate in and out of the control volume. Exergy 
destruction is expressed in terms of product and fuel for a 


specific component. 
Ex = Èr k m Epr-Enn 


(3) 


The exergy efficiency,w,, and the exergy destruction ratio is 
equally defined for the ktt component as: 


(4) 
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Furthermore, Eqs (1) to (5) were applied to Figure 1 to derive 
the energy quantities for the system presented in Table 1 in 
Appendix. 


3.1.1 Exergoeconomic modeling 
The components’ purchase and equipment cost (PEC) 
are written as functions of their operating parameters. The 


general cost balance for a control volume for the kt” 
component is presented [21]. 
Can + LiCin + Ze = Delon t wk (6) 


Where Coa, is the cost associated with the j^" sum of exergy 
streams to the system's k^, Y; Ci, is the levelized cost rate 
for the kt” component, Z,. The cost associated with the jt” 
sum of exergy streams from the system's kt” component 
is Ye Caki while the work associated with the kt” component 
in the system is denoted with the term C uk: 

The cost of the jt? stream is related to the cost of specific cost 
and the exergy, work, or heat with the relationships: 
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Cw k = CyW (8) 
Cox = qQ (9) 


The cost rate Z;, for the components is present as: 


> _ PECFXCRFX$ 

Ze = — NX3600 (10) 

PECF, CRF, and ¢ represent the purchase of equipment cost 

function, capital recovery factor, and maintenance factor, 

respectively. The system's annual operational hours are 

denoted with N while the capital recovery factor is expressed: 
ifati|™ 


CRF = aaa (11) 


Where n is the estimated plant life in years, another exergy- 
related index is the cost of exergy of product and fuel for the 
component. A detailed description of the component cost of 
fuel and product is found in [22]. The general relationship for 
evaluating the specific cost of a product for the kth 
component cp, ($/kJ), and that of fuel cp x ($/kJ), as well as 
the exergoeconomic factor, f is presented in Eqs (12) - (14), 
respectively. 


Cp ($/k]) = zt (12) 
crx ($/K) = F (13) 
f= bum J 


The cost of exergy destruction is expressed in Eq (15). A 
summary of the cost-related terms for Figure 1 is presented 
in Table 2, while the PEC for the plant components is shown 
in Table 3. 


Cp, = Cp kEp,k (15) 


3.2 Thermoenvironmental analysis 

The environmental effect is evaluated by approximating 
the measure of pollutants produced by the plant. These 
include nitrogen oxide, minox (kg/s), CO2, mco, (kg/s) and CO, 
mco (kg/s). The quantity of these emissions produced and 
their production rates depend on the following indices: 
retention time, T (s), combustion chamber pressure drop, 
APCC (kPa), adiabatic flame temperature, Tj; presented [23, 
24]. Consequently, the emission rates of production and the 
harmful emission factorFzr, are defined in Eqs (16) to (20). 


1.5x 10155056 - (7110/Tpz) 


Two, ges (aPeoy (16) 
: 1.79x109 7059 (7800/Tpz) 

Mco p, ^os, (22E0) (17) 
. = CO . Mcoz 18 
Mco, = YOUZMg mg (18) 
Coz sp = 3600 (G2) (19) 
E = Mnox*Mco*Mco, (20) 


mg 


Where COz;,(kgco,jmwn) describes the amount and the 
specific CO; emissions. Similarly, m,, M , yCOz and Mco, are 
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the rate of mass flow of flue gas, flue gas molar mass, mass 
fraction and the CO, molar mass. The adiabatic flame 
temperature Tpz is defined by Eq (21). All parameters, 
constants and the terms x, y and z in Eq (21) are estimated 
according to the procedure in [24]. 


Tpz = Ac*exp[ p(o + A)? |n* 0Y y? (21) 


3.3 Sustainability indicator 
3.3.1 The exergetic utility index 

The exergetic utility index (EUI) is depicted in Eq (22) 
and measures the extent of exergy resource utilization in a 
system regarding the same system's net output. Itis a function 
of combustion efficiency,Acc, net output, Wnet, exergy input, 
Ein and the exergy efflux to the environment, É,,,, [15]. 


EUI ettet. (22) 
Exin—Exout 
3.3.2 Exergo-thermal index 
The exergo-thermal index (ETI) measures the thermal 
impact of the system on the environment during the energy 
conversion process. Low values of ETI are desired and can be 
achieved by continuously utilizing high-temperature flue gas 
from energy conversion systems in powering other low-heat 
bottoming cycles, expressed as: 


I AccXWnet x1 
(Exin Éxout) Y 


ETI = ES (23) 
Where Y denotes the enviro-thermal conservation factor, 
expressed as: 


=~ 
te (24) 
Where Ty is the ambient temperature and T; is the 
temperature of the flue gas. 


3.3.3 Exergetic sustainability index 

The exergetic sustainability index (ESI) compares the 
magnitudes of a system’s net product to its exergy 
destruction. A system with ESI less than unity underutilizes 
the exergy resources, and it is not desirable. Conversely, a 
system with ESI greater than unity is desirable as its net 
output exceeds the total exergy destruction. The ESI is 
expressed by [25]. 


Wnet 


ESI = 4 (25) 
ED Total 

4. Results and discussion 

4.1 Thermodynamic properties and operating 


conditions 

The simulation results of the system were obtained using 
a developed computer code written in EES. The initial 
operating conditions are presented in Table 4. The model 
development was built on the succeeding assumptions: the 
surrounding temperature and pressure exist at 25 °C and 
1.013 bar, respectively. The system and its subcomponents 
were evaluated at steady-state conditions. Pressure and 
temperature variations were neglected. The system's 
boundaries were treated as adiabatic. The working fluid for 
the Kalina system is ammonium water solution at 0.28. 
Likewise, the properties of the system, including the enthalpy 
and exergy flow rates used to calculate system performance, 
are shown in Table 5. 
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Table 2. Component cost and auxiliary equations for the multigeneration plant 


Component Exergoeconomic balance Auxiliary equation 
GT LPC C, wpe t Zupc = C2 Nil 
GT intercooler Cot Cag + Zinr = Cz + Cao Ġ È — CE, = 0 
GT HPC C3 + Cic + Zypc = Cy Nil 
GT Hex Cyt Cat Zugy = Čs t+ Cp CÉ;- CE, = 0 
GT.CÉ (bU E =Ġ Nil 
GT HPT C; Z uper = Ca Cwupe + Owe Ès — GE; =0 
Cy, pc Wape — Curgpc Wipc = 0 
GT REH Cy Ate Tus = Ġo Nil 
GT LPT Ġo + Zipp = Cr+ Ĉwipr CroF a -= Cii Eio -0 
Kal. Vap Gen. [um + Caz + Zyg = = [um Tg CE — Ci3E12 = 0 
Kal. Sep 1 Gs + Loup, = = Gig t C5 C;gE5; — C22E18 = 0 
Kal. Valve 1 Cj; + Zy, = Cog Nil 
Kal. Valve 2 Cj + Zy, = C Nil 
Kal. Turb. Cis + Zkatrury = Coo Cw py CWeump,  ĈisË20 — CaoFig = 0 
Turb s 
Wratrurb Wpump, — Cwpump, W aig =0 

Kal. Sep 2 Coo Zsgp, = C21 + Cog Ca, Eg — CogE5; = 0 
Kal. CND 1 Cog + Cso + Zenpi = C29 + Coo C28E29 — C29E28 = 0 
Kal. Valve 3 Coo Zy, = Cao Nil 
Kal. Evap 1 Czy + Cop + Zgyp = Car + Cog Cs; Egg — CogEe7 = 0 
Kal. HEX 1 Cog + C36 + Zugx = 62; + C37 Cy6E27 — C27E26 = 0 
Kal. HEX 2 Cs; + Cas + Zuex = [om + Cae C32E33 — C33E32 = 0 
Kal. CND 2 Cas + Cs; + Lenn, = Čs4 + Čss C33E 34 — C34E33 = 0 
Kal. Pump 1 C34 + CW pump, + Zpump, = Cas Nil 
VAS desorber Cos + Cag + Zpgsg = Cog + Cag + Cas CosE 26 — CopE 2s = 0 

. i . . CagEas = Cas E 3g =0 
VAS HEX 3 C43 + Cast Zoey = Cag + Cae C45 E4g — C4gE45 = 0 
VAS Valve 4 Cas + Zy, = Cao Nil 
VAS Pump 2 Cy + en; + Zpump, = Cas Nil 
VAS Absorber Car + Cay + Cos + Zapsg = Car + Coe Css — Cos| * £j; — Caz * [Eso — Es] = 0 
VAS EVP 2 Cao + D + Zeyp = Cu + Cog Cu Eg, — Cu Es = 0 
VAS Valve 5 Cag Zy, = Caz Nil 
VAS CND 3 Coq + n + Zens = = Ct Cs Caglzg — CagEgg = 0 
Water heater Cc Cis + Zwy = = Cyt Cre CE, — CES = 0 


PEM Electrolyser C5 + Cig + Zevecrr. = Cs3 + Csa 


PEM Hex Coat Coot Coo + Zugy = Csi + Coo 


PEM Ozseparator C54 + Zoxygensey = Css + Css 


[nm + Cs ||Éss + Esa|— ICs3 * Css | Es + Ey9| =0 
Css. + Cool, * Ego — Cog * [Ese + Ecol =0 
CesEse — Cs Egg = 0 


4.2 Performance analysis of the system based on 

operating data 

The performance indices considered for the study are 
presented in Table 6. The system net energy and exergy 
efficiency was calculated at 53.48 and 50.05 %, respectively. 
The multigeneration system improved the topping cycle by 
19.03 and 2.58 % in energy and exergy efficiencies, 
respectively. The improvement was ascribed to the added 
products from the bottoming cycles estimated at 30,178 kW. 
These additional products include the power generated from 
the Kalina turbine, VAS and Kalina system cooling, and the 
quantity of hot water produced from the domestic water 
heater. The net power from the topping cycle stood at 55.605 
MW. Also, the system hydrogen production rate stood at 
0.1524 kg/h, with about 3498 kW of cooling and a coefficient 
of performance of 4.304. 


4.3 Component exergy destruction (ED) rates 

The system components’ ED rates are shown in Table 7. 
About 85.81 96 of EDs were from the GT topping cycle due to 
significant ED rates around the combustion chamber (CC). 
The GT CC and reheater contributed approximately 54.65 and 
22.57 96 of the ED, respectively, to the total ED. The main 
reason for the large ED in the CC is the large temperature 


difference between the combustion gases and the hot air. 
Using alternative preheated fuels before combustion will 
reduce ED in the CC. In the present study, operating the 
reheater with biomass syngas at 154 °C reduced ED by nearly 
25 96. Furthermore, with the integration of the bottoming 
cycles, about 15089.77 kW ED was avoided, equivalent to 
14.19 96 of the total ED. Also, an additional product of 30178 
kW was achieved, which increased the  exergetic 
sustainability of the bottoming cycle to approximately 1.99. 
The latter validates retrofitting the topping cycle, the Kalina 
and VAS cycles. 


4.4 Exergoeconomic parameters of the plant 

The values of the initial investment, monetary flow rate 
and levelized capital cost (LECC) rate are shown in Table 8. 
The GT HPT hasthe highest levelized cost per hour, calculated 
at 41.5 $/hr, followed by GT LPT and GT LPC, estimated at 
34.19 and 12.17 $/hr, respectively. Similarly the 
exergoeconomic parameters of the subsystems are depicted 
in Tables 9-12. From Table 9, the combustion chamber and 
reheater contribute about 39.903 and 20.613 $/hr to the total 
ED cost calculated at 106.4 $.hr. The CC and reheater have the 
least exergoeconomic factors of 0.112 and 0.279, respectively. 
One method to reduce ED cost in the CC is to improve 
combustion efficiency and reduce the temperature difference 
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between the working fluid and combustion fuel. However, 
this can be achieved by firing the CC with preheated fuel. In 
the Kalina subsystem (Table 10), the vapour generator (VG) 
cost improvement potential is primarily dependent on the 
equipment cost rates than the cost due to ED. 
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Subsequently, since the VG cost is directly linked to the 
quantity of heat transfer, the choice of working fluid and 
efficient design of the heat transfer area is germane in 
reducing the cost of the system components. 


Table 3. Component cost functions, cost of product and cost of fuel [15, 18, 22] 


Component Cost function Cost of product Cost of fuel 
I$] 
GT LPC 395m, [Pal |P 6-6 X Wupr 
0.9 — n, pc IP, P, . : . : 
GT intercooler ae Ayr [279 C49 — Cag C; — C 
0.093 
GT HPC 3955, [Pal [Pa C=C [1 = x|. Wupr 
0.9 — nupc |P3 P 
GT Hex 138 Ayuxx [279 6 — Cy Cu - C 
0.093 
GT CC 6 Cs + Cs 
46.081; 
P |1 + exp(0.018T; — 26.4)| 
0.995 — P 
5 
GT HPT 479.34m; P, W, 6, — 6 
In |<] |1 + exp[0.036T; — 54.4]] DET NL 
0.92 — quer Pg à 7 
GT REH [o Cs t C 
46.08rhg 2 diiit: 
|1 + exp(0.018T,, — 26.4)| 
0.995 — Pro 
Pg 
479.34m P. ó 
elie 7 | in [29] |1 + exp[0.0367,, — 54.4] | Wier Cao Cn 
0.92 — Mer 11 . 
Kal. Vap Gen. 33m AvAscEN C, — C37 C; — Chg 
0.093 
Kal. Valve 1 37 Pj s C23 C22 
P23 
Kal. Valve 2 an Py pee on Coy 
Poa 
479.34m P i | j mia C 
Kal: Turb. ee! m [21 4 + exp[0.036T,s — 54.4] | Mier + Wratp,  Wyasy, Cis = Coo 
0.92 — "nr Po . : . . 
Kal. CND 1 — Qxatcnp1 Cso — Cso Cog — Cs 
` 0.15AT Kaicnp1 
Kal. Valve 3 T Pro pee C30 Cs 
P30 : 
Kal. Evap 1 X Qvoisvsi o8» C31 — Cao C5; — Ces 
0.15ATkaigyps 
Kal. HEX 1 130 [Arten Cs; — C36 Cas — Cas 
0.093 
Kal. CND 2 ere Üaicxpa Csg — Cs; C33 — C34 
0.15ATkaicnn2 ] . 
Kal. Pump 1 705.5 [0.061Wca 071 0.2 Css — C34 Wratp, 
Pump1 1 m Np 
VAS HEX 3 130 (Av4snexs 5 Cas — Cas Cas — Cae 
0.093 
VAS Valve 4 37 P» idi Ca; Cas 
40 
i 0.71 0.2 = j 
VAS Pampa 705.5 |0.001Wratpumps Pico Ton 
p 1-1» 
A 0.85 ^ 
VAS EVP 2 T 4l Opisz Cu — Cao Cos — Coa 
` 0.15ATy4s EVP2 
VAS CND 3 516.62 QvAscND3 Cs; — Cai Csg — C39 
0.15ATy4scu pa E 
Water heater 130 |Attaren Cig — Cis C; — Cu 
0.093 


45 


FI Abam et al. /Future Technology 


Table 4. Design operating data [15, 23] 
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Parameter Unit Value 
Ambient temperature oC 25 
Ambient pressure Bar 1.013 
GT Lower compression ratio Dim. 3.162 
GT Higher compression ratio Dim. 3.162 
Overall pressure ratio Dim. 10 
GT Heat exchanger effectiveness % 75 
Low pressure turbine isentropic efficiency % 85 
High pressure turbine isentropic efficiency % 85 
Low pressure compressor isentropic efficiency % 80 
High pressure compressor isentropic efficiency % 80 
mass of air to the topping cycle kg/s 200 
mass of gas at the inlet to the first combustion chamber kg/s 3.131 
mass of gas at the inlet to the reheater kg/s 21.49 
High pressure turbine inlet temperature oK 1350 
Low pressure turbine inlet temperature oK 1200 
The exit temperature of intercooler water oC 85 


The pump has an exergoeconomic factor of 1, suggesting a 
zero contribution to ED costs. This is realistic as the system 
pumps operated in isentropic conditions. The total cost of 
product and fuel for the GT (Table 9) is calculated at 2.67 and 
6.06 $/GJ, respectively, while that for the Kalina (Table 10) 
and VAS subsystems (Table 11) exist at 32.92 and 84.89 $/GJ, 
and 17.97 and 25.93 $/GJ, respectively. The results show a 
lower average energy cost from the GT (0.836 $/GJ) 
compared to the Kalina subsystem and the VAS. For the PEM 
electrolyzer, economic parameters are equally presented in 
Table 12. Furthermore, the exergoeconomic factors (fk) for 
all the subsystems are similarly presented in Tables 9-12, 
estimated at 59.59%, 68 96, 8.656%, 30.71%, 47.41% for GT, 
KAL, VAS, and PEM electrolyzer respectively. Low values of 
fk for any system component indicate high ED cost, which 
signifies a high prospective for improvement, while high 
values of relative cost difference (rk) indicate prospects for 
system optimisation. The rk for the GT, KAL, VAS, and PEM 
electrolyzer subsystems are calculated at 127 96, 157 96, 30% 
and 30.71 96, respectively. The joined effects of rk and fk 
show that the GT and KAL subsystems have the highest 
potential for improvement, followed by VAS and PEM 
electrolyzer, which have the narrow potential for 
optimisation. 


4.5 Sensitivity analysis 
4.5.1 Effect of dead state temperature on system 

efficiencies 

Figure 2 presents the effect of dead state temperature 
on system performance. The energy and exergy efficiencies of 
the GT and Kalina cycles and the coefficient of performance of 
the VAS were observed. The dead state temperature (DST) 
ranged between 288 and 302 K. The results show that the 
DST increase led to a reduction in system efficiency. At high 
DSTs, the compressor work increases, reducing the net 
system output, especially at the topping GT cycle. The results 
also show that for every 5 Krise in DST, the energy and exergy 
efficiency reduced by 0.40 and 0.44 96, respectively. The COP 
of the VAS was estimated at 4.30 and is constant throughout 
the DST range, showing that the DST has no direct link with 
COP measured parameters. 


4.5.2 Effect of heat input on the PEM electrolyzer heat 

exchanger 

The effect of the required temperature of heating at 
the inlet to the PEM heat exchanger (HEX) is investigated on 
the hydrogen yield rate (Figure 3). The water inlet from the 
condenser of the Kalina system improves the quantity of the 
hydrogen production rate due to higher PEM water inlet rate 
to cater for the high energy exchange required in the HEX. 
Since the temperature of the water entering the PEM 
electrolyzer is fixed at 95°C, the figure demonstrates the 
range of water temperatures at the PEM HEX. At most, 
102.02°C is required for optimum hydrogen yield with other 
operating PEM electrolyzer parameters constant as indicated 
in its modelling. 


4.5.3 Effect of gasifier temperature on hot water rate 

power output 

Figure 4 shows the effect of gasifier temperature in the 
topping cycle on the quantity of hot water production at a 
given hot water temperature and the magnitude of power in 
the Kalina turbine. Very high gasifier temperatures result in 
slightly reduced biomass syngas calorific value. Thus, the 
energy level of the bottoming cycles, including the water 
heater, is reduced, affecting the quantity of hot water at 959C. 
However, the power output from the Kalina turbine was 
constant at 181 kW. This is attributable to the fixed operating 
conditions of the Kalina system within a range of heat transfer 
requirements around the Kalina vapour generator. 


4.5.4 Effect of ammonia mass fraction on cooling and 

turbine output 

The effect of ammonia mass fraction on the power 
output and cooling rate of the Kalina system is presented in 
Figure 5. In the design condition, a mass fraction of 0.28. The 
result indicates that increasing ammonium mass fraction 
leads to a higher turbine work and cooling rate. However, at 
mass fractions higher than 30 per cent, and with the operating 
pressures in the system, the pump work requirements are 
significantly higher than the turbine output. Therefore, it is 
not feasible to operate the system with an ammonia mass 
fraction in excess of 0.30. 
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Table 5. Thermodynamic properties at the state points 


State Temperature Pressure Enthalpy Entropy Mass Exergy 
[°c] [Bar] [kJ/kg] [kJ/kg-K] [kg] [kw] 
1 25 1.013 298.40 5.695 200 0.00 
2 170.1 3.203 445.00 5.765 200 25150.00 
3 120 3.203 397.90 5.643 200 22965.00 
4 311.3 10.13 590.90 5.72 200 57016.00 
5 523.2 10.13 818.10 6.051 200 82706.00 
6 25 10.13 42300.00 - 2.533 124996.00 
7 977 10.13 1337.00 6.565 202.5 157799.00 
8 689.7 2.843 1004.00 6.628 202.5 86580.00 
9 427 2.843 19841.00 - 2.991 62661.00 
10 927 2.843 1278.00 6.882 205.5 128628.00 
11 696.6 1.013 1012.00 6.932 205.5 70840.00 
12 484.7 1.013 776.00 6.658 205.5 39178.00 
13 471 1.013 761.10 6.638 205.5 37326.00 
14 349.6 1.013 631.20 6.448 205.5 22305 
15 20 1.013 83.30 0.294 85 0.00 
16 95 1.013 397.40 1.248 85 2523.00 
17 180 20 1264.00 3.639 2.5 862.60 
18 180 20 2164.00 5.901 0.97 606.80 
19 - - - - - 4.774 
20 136.9 7 1971.00 5.901 0.97 419.60 
21 136.9 7 509.50 1.753 0.0644 4.69 
22 180 20 694.10 2.205 1.53 206.70 
23 79.8 0.8 694.10 2.39 1.53 122.20 
24 77.6 0.8 509.50 1.821 0.0644 3.39 
25 79.7 0.8 686.70 2.367 1.594 125.60 
26 60 0.8 176.80 0.8905 1.594 14.33 
27 55 0.8 146.50 0.799 1.594 9.42 
28 136.9 7 2080.00 6.209 0.9056 412.80 
29 42.9 7 -42.60 0.48 0.9056 36.90 
30 -5.5 0.8 -42.60 0.5609 0.9056 15.08 
31 15 0.8 332.40 1.909 0.9056 -9.22 
32 43.3 0.8 213.80 1.236 2.5 25.78 
33 38 0.8 123.00 0.9479 2.5 13.61 
34 25 0.8 -73.38 0.3073 2.5 -0.06 
35 25.1 20 -71.24 0.3073 2.5 5.28 
36 40 20 -7.27 0.5168 2.5 9.19 
37 50 20 35.55 0.6514 2.5 15.91 
38 79.7 0.07424 2649.00 8.478 1.5 189.40 
39 41 0.07424 171.10 0.5836 1.5 2.20 
40 1.7 0.006812 171.10 0.6229 1.5 -15.35 
41 1.7 0.006812 2503.00 9.114 1.5 -312.80 
42 34.6 0.006812 92.00 0.1987 11.16 3.02 
43 34.6 0.006812 92.00 0.1987 11.16 3.02 
44 67 0.006812 156.30 0.3978 11.16 58.69 
45 79.7 0.006812 229.30 0.4167 9.661 40.51 
46 45 0.006812 169.30 0.2361 9.661 -19.08 
47 35 0.006812 169.30 0.1817 9.661 137.40 
48 20 1.013 83.30 0.294 37.58 0.00 
49 80 1.013 334.30 1.073 37.58 705.30 
50 20 1.013 83.30 0.294 5 0.00 
51 95 1.013 397.40 1.248 5 148.40 
52 85 1.013 355.30 1.132 10 221.60 
53 85 1.013 4843.00 56.17 0.0000423 0.00 
54 85 1.013 355.30 1.132 5 110.80 
55 85 1.013 55.01 0.1683 0.000336 0.00 
56 85 1.013 355.30 1.132 5 110.80 
57 20 1.013 83.30 0.294 8.383 0.00 
58 34 1.013 141.90 0.4893 8.383 3.18 
59 20 1.013 83.30 0.294 7.065 0.00 
60 85 1.013 355.30 1.132 7.065 156.60 
61 20 1.013 83.30 0.294 44.42 0.00 
62 40 1.013 166.96 0.5702 44.42 59.99 
63 25 1.013 298.42 5.695 151.5 0.00 
64 2 1.013 275.30 5.614 151.5 142.40 
65 20 1.013 83.30 0.294 69.55 0.00 
66 35 1.013 146.00 0.5029 69.55 35.41 
67 25 1.013 298.40 5.695 13.01 0.00 
68 -1 1.013 272.30 5.603 13.01 15.73 
69 33 1.013 137.80 0.4762 15.45 4.14 
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Table 6. Performance parameters of the energy system 


Performance Index Unit Value 
LPC Power MW 29.151 
HPC Power MW 38.743 
LPT Power MW 55.605 
HPT Power MW 68.754 
GT Net Power Output MW 55.605 
GT Relative Exergy Efficiency 96 29.63 
Kalina Relative Exergy Efficiency 96 11.13 
Kalina Relative Energy Efficiency 96 16.97 
VAS COP - 4.304 
Kalina turbine power kW 187.2 
Kalina pump power kW 5.337 
Kalina evaporator cooling kW 339.6 
Kalina Exergy of cooling kW 24.3 
Kalina net power kW 181.9 
VAS evaporator cooling kW 3498 
Exergy of cooling kW 297.4 
VAS desorber heat kW 812.9 
PEM hydrogen output kg/h 0.1524 
System net output MW 86.320 
System net input MW 161.403 
System net energy efficiency 96 53.48 
System net exergy efficiency 96 50.05 


4.5.5 Gasifier mass flow rate on thermo environmental 

parameters and hot water production rate 

The effect of the gasifier mass flow rate at the reheater 
was investigated on the quantity of hot water production and 
thermo-environmental parameters, as shown in Figure 6. The 
result indicates that the gasifier mass flow rate increases the 
amount of hot water production when maintained at 95 °C. 
This slight increase in gasifier mass flow rate has a negligible 
impact on ETI, decreasing slightly from 2.115 to 2.109 units. 
The trend is attributed to the variation in the flue gas 
temperature after the water heater. Also, similar 
thermodynamic conditions are responsible for a low 
comparative variation on the EUI when the gasifier mass flow 
rate increases. Conversely, increasing the gasifier mass flow 
rate reduces the overall ESI. The system's ESI dropped from 
0.822 to 0.358, which is attributed to large exergy destruction 
in the reheater and LPT. 


4.5.6 Effect of ammonia mass fraction on thermo- 

environmental parameters 

Figure 7 shows the effect of ammonia mass fraction on 
thermo-environmental parameters. The variations in the 
ammonia mass fraction slightly affect the ETI, EUI, and ESI. 
Higher ammonia concentration results in a rich and increased 
quantity of ammonium vapour for expansion in the turbine 
and evaporation in the evaporators. Thus, the power output 
ofthe Kalina system rises. Nonetheless, since the EUI, ETI, and 
ESI are wholly affected by the variation of the net output from 
the entire multigeneration plant, the effect of Kalina output is 
negligible on these thermo-environmental parameters. 
However, due to high exergy destruction in the Kalina vapour 
generator resulting from the increase in the ammonium 
concentration, the ESI reduces slightly from 0.5279 to 0.5251 
between ammonia concentrations of 0.25 and 0.40. 


4.5.7 Effect of primary zone temperature on 
environmental emissions 
Figure 8 present the effect of the primary zone 
temperature on the production ofCO, CO? and NOx 
emissions. The parameters of the primary zone temperature 
in the combustion chamber were estimated using the 
methodology in [24]. 
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Table 7. Summary of systems component and total exergy 
destruction 


Component Exergy [96] [96] 
destruction ED/cycle total of 
(ED) and ED 
[kW] component 
GT CC 49903.00 54.652 46.901 
GT HEX 5971.00 6.539 5.612 
GT HPC 4692.00 5.139 4.409 
GT HPT 2466.00 2.701 2.318 
GT INTCL 1481.00 1.621 1.392 
GT LPC 4001.00 4.382 3.760 
GT LPT 2183.00 2.391 2.052 
GT reheater 20613.00 22.575 19.373 
GT total ED 91310.00 100.000 85.818 
KAL. CND1 219.3.00 15.517 0.206 
KAL. CND2 10.480 0.742 0.009 
KAL. EVP1 8.568 0.606 0.008 
KAL. HEX1 1.813 0.128 0.002 
KAL. HEX2 8.265 0.585 0.008 
KAL. PUMP1 0.003 0.00018 2.5E-06 
KAL. SEP1 49.040 3.469 0.046 
KAL. SEP2 2.147 0.151 0.002 
KAL. TURB 0.000 0.000 0.000 
KAL. VALVE1 84.560 5.983 0.079 
KAL. VALVE2 1.301 0.092 0.001 
KAL. VALVE3 21.820 1.544 0.021 
KAL. VG 1006.00 71.181 0.945 
Kalina total ED 1413.290 100.00 1.328 
VAS absorber 213.900 21.543 0.201 
VAS condenser 2 127.200 12.811 0.119 
VAS desorber 318.800 32.108 0.299 
VAS evaporator 2 155.100 15.621 0.146 
VAS HEX 3 3.942 0.397 0.004 
VAS pump 2 0.000 0.000 0.000 
VAS valve 4 17.550 1.768 0.016 
VAS valve 5 156.400 15.752 0.147 
VAS ED 992.890 100.000 0.933 
Water heater 12502.00 100.000 11.750 
PEM electrolyzer 115.600 63.658 0.109 
PEM HEX 65.990 36.338 0.062 
PEM oxygen 0.0060 0.0032 5.5E-06 
separator 


The primary zone temperature affects the emission rates 
at critical temperatures in excess of 1895 °K. The rate of NOx 
emissions are higher (10.57 kg at 2400 °K) compared to CO 
(0.01380 kg at 2400 »K). Consequently, it is recommended to 
operate the system at optimum primary zone temperature to 
reduce NOy emissions. The CO2 emissions from the system 
are not directly related to the primary zone temperature but 
to the emissions coefficient of the natural gas after 
combustion and its mass flow rate. However, this system 
calculated the instantaneous CO? emissions at 27.35 kg/s ata 
flue gas mass flow rate of 205.5 kg/s. 
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Table 8. Summary of initial investment, monetary flow rate 
and Levelized capital cost rate 


Levelized 
Plant component Purchase Levelized cost per 
of cost per hour 
equipment year ($/hr.) * 
cost ($) * ($/yr.) 
GT CC 213978 38022 5.038 
GT HEX 155098 27559 3.652 
GT HPC 517708 91992 12.19 
GT HPT 1762604 313198 41.5 
GT INTCL 93214 16563 2.195 
GT LPC 517708 91992 12.19 
GT LPT 1452180 258038 34.19 
GT reheater 213093 37864 5.017 
KAL. CND1 1953 347 0.04598 
KAL. CND2 1953 347 0.04598 
KAL. EVP1 1953 347 0.04598 
KAL. HEX1 3353 595.9 0.07895 
KAL. HEX2 3353 595.9 0.07895 
KAL. PUMP1 5405 960.5 0.1273 
KAL. SEP1 13067 2322 0.3076 
KAL. SEP2 13067 2322 0.3076 
KAL. TURB 156801 27862 3.692 
KAL. VALVE1 325.5 57.83 0.007663 
KAL. VALVE2 325.5 57.83 0.007663 
KAL. VALVE3 325.5 57.83 0.007663 
KAL. VG 191786 34079 4.515 
VAS absorber 1060 188.3 0.02495 
VAS condenser 3 1137 202.1 0.02678 
VAS desorber 2679 476.1 0.06308 
VAS evaporator 2 1017 180.7 0.02395 
VAS HEX 3 1254 222.8 0.02952 
VAS pump 2 5321 940.3 0.0123 
VAS valve 4 325.5 57.83 0.007663 
VAS valve 5 325.5 57.83 0.007663 
Water heater 2444 434.3 0.05755 
PEM electrolyser 3823 679.2 0.0225 
PEM HEX 955.6 169.8 0.0225 
PEM oxygenseparator 1147 203.8 0.027 


4.5.8 Effectofgasifier temperature and mass flow on the 

unit cost of electricity 

The effect of gasifier temperature and mass flow rate 
on the unit cost of electricity is shown in Figure 9. The 
temperature of the gasifier and the mass flow correlate with 
theunit cost of electricity (UCOE) regarding the system power 
output. The UCOE based on the initial design conditions was 
obtained as 21.44 N/kWh (USD 0.0329/kWh), which is more 
attractive than the current tariff plan in Nigeria of 54 N/kWh 
(USD 0.083/kWh). The gasifier temperature and the mass 
flow rate increase leads to a decrease in the UCOE following a 
high output generation from the system. Therefore, the 
reduction in UCOE is more with respect to the gasifier 
temperature increase than the gasifier mass flow rate. 


4.5.9 Effect of ammonium mass fraction on the unit cost 
of electricity 
The effect of ammonia mass fraction on the UCOE and 
PEC for the entire system is shown in Figure 10. The increase 
in ammonia mass fraction leads to a slight increase in the 


Hot water production rate (kg/s) 
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UCOE. High ammonium mass fraction results in high turbine 
output, which increases the turbine purchase equipment cost. 
Although there is an additional power generation from a high 
ammonium mass fraction, the UCOE increase is negligible. 


—a— Energy efficiency 
—4—VAS COP 


—o— Exergy efficiency 
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Figure 2. Dead 
parameters 
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Figure 3. Effect of heat input on the PEM electrolyzer heat 
exchanger 
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Figure 4. Gasifier temperature effect on hot water production 
rate and power output 
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Table 9. Exergoeconomic parameters for the GT and reheater topping cycle 


Component (E; Ga Ep Cp Z Z+Cp f Tk 
($/G)) ($/GJ) (MW) ($/hr.) ($/hr.) ($/hr.) (%) (%) 
GT CC 0.2217 0.3007 49.903 39.8243 5.038 44.8623 11.30 35.63 
GT HEX 0.2139 0.3028 5.971 4.597719 3.652 8.2497 44.30 41.55 
GT HPC 0.3649 0.6516 4.692 6.165175 12.19 18.355 66.40 78.53 
GT HPT 0.3007 0.4851 2.466 2.669287 41.5 44.169 93.40 61.34 
GT INTC 0.6967 3.0239 1.481 3.714362 2.195 5.909 37.10 334.06 
GT LPC 0.4851 0.6969 4.001 6.986962 12.19 19.177 63.50 43.67 
GT LPT 0.2139 0.3930 2.183 1.680655 34.19 35.871 95.30 83.77 
GT REH 0.1749 0.2139 20.613 12.98443 5.017 18.001 27.90 22.22 
Table 10. Exergoeconomic parameters for the Kalina power-cooling bottoming cycle 
Component Cr Cp Ep (55 Z Z+Cp fg Fi 
($/GJ) ($/G)) (MW) ($/hr.) ($/hr.) ($/hr.) (%) (%) 
KAL. CND1 2.5492 6.1994 0.2193 2.0126 0.0460 2.0585 2.230 143.190 
KAL. CND2 2.4799 14.579 0.0105 0.0935 0.0459 0.1395 32.95 487.876 
KAL. EVP1 1.0000 2.1805 0.0086 0.0000 0.0459 0.0459 100.00 118.051 
KAL. HEX1 1.9857 4.7164 0.0018 0.0130 0.0789 0.0919 85.89 137.515 
KAL. HEX2 2.4787 13.327 0.0083 0.0738 0.0789 0.1527 51.70 437.671 
KAL. PUM1 8.0153 14.968 2.6E-06 7.5E-05 0.1273 0.1274 99.94 86.738 
KAL. SEP1 2.1012 2.3329 0.0490 0.3710 0.3076 0.6786 45.33 11.025 
KAL. SEP2 2.3329 2.5496 0.0022 0.0180 0.3076 0.3256 94.46 9.2890 
KAL. TURB 2.3326 7.5825 0.0000 0.0000 3.6920 3.6920 100.00 225.064 
KAL. VAL1 2.3329 3.9644 0.0846 0.7102 0.0077 0.7179 1.07 69.928 
KAL. VAL2 2.5474 4.1527 0.0013 0.0119 0.0077 0.0196 3.91 63.019 
KAL. VAL3 2.5542 6.3918 0.0218 0.2006 0.0077 0.2083 3.68 150.247 
KAL. VG 0.2129 1.9493 1.0060 0.7713 4.5150 5.2864 85.41 815.217 
Table 11. Exergoeconomic parameters for the vapor absorption bottoming cycle 
Component Cr C Ep Cp Z Z+Ċp fx T 
($/GJ) ($/GJ) (MW) (S/hr.) ($/hr.) ($/hr.) (90) (%) 
VAS ABS 0.0909 1.2461 0.2139 0.0701 0.0250 0.0950 26.26 1269.62 
VAS CON 3 1.4121 4.5304 0.1272 0.6466 0.0268 0.6734 3.980 220.82 
VAS DESB 0.8303 1.4123 0.3188 0.9529 0.0631 1.0160 6.210 70.097 
VAS evaporator 2 1.0001 1.0224 0.1551 0.0000 0.0240 0.0240 100.00 2.23 
VAS HEX 3 0.3313 0.6916 0.0039 0.0047 0.0296 0.0342 86.26 108.76 
VAS pump 2 15.014 17.473 0.0000 0.0000 1.02E-07 1.02E-07 100.00 16.38 
VAS valve 4 0.1408 0.3406 0.0176 0.0089 0.0077 0.01660 46.28 141.91 
VAS valve 5 20.663 24.9180 0.0031 0.2290 0.0050 0.2340 97.69 20.59 
Table 12. Exergoeconomic parameters for the PEM electrolyzer 
Component Gr Cp Ey Cp Z Z+Cp fx Tk 
($/G]) ($/G]) (MW (S/hr.) ($/hr.) ($/hr.) (96) (96) 
WATER heater 0.2140 1.2793 12.502 9.6297 0.05755 9.6873 5.9 497.93 
PEM Electrolyzer 4.5650 9.3813 0.116 1.8998 0.0225 1.9223 1.17 105.50 
PEM HEX 6.3660 6.7198 0.066 1.5124 0.0225 1.5348 1.47 5.56 
PEM oxygen separator 9.3813 9.4491 5.79E-06 0.0002 0.0270 0.027196 99.28 0.723 


4.5.10 Effect of gasifier temperature and mass flow on 

system total cost 

In Figure 11, the effect of gasifier mass flow rate and 
temperature is investigated on the system's total cost. The 
system's net cost was estimated at $ 5335411.5 based on the 
operating conditions in Table 6. The total cost comprises 
purchasing equipment (PEC) for all system components. Thus 
an increase in both the gasifier temperature and mass flow 
rate significantly raises the total PEC of the system. 


The resultant effect of increased gasifier temperature and 
mass flow rate cumulatively enhances the power 
requirements of the LPT. Consequently, the PEC of the LPT, 
GT HEX, Kalina vapour generator, and water heater 
contributes about 33.76 96 of the entire plant's PEC. 


4.6 Optimum parameters 
For optimal performance, the study considered the 
following objective functions (OBF): total exergy efficiency 
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(V) and the total product cost rate. The yy is to be maximized, 
while the total cost rate is to be reduced or maximized. The 
total cost rate model equation of the multigeneration plant is 
integrated with the cost rate of the pollution damage. 
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Figure 5. Effect of ammonia mass fraction on Kalina cycle 
performance 
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Figure 6. Gasifier mass flow rate on environmental 
parameters and hot water production rate 
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Figure 7. Effect of ammonia mass fraction on thermo 
environmental parameter 


The OBFs are described as follows: 


Woverail = 

à i . . . fg : 6 jp 

[Ew or Ew prt Éwrier Eneas ÉRvas {EwrectEwupct Zier Ew] 
Aa 


(26) 


UCOE (N/kWh) 


2 


UCOE (N/kWh) at diff. gasifier temp. 


February 2024| Volume 03 | Issue 01 | Pages 40-55 


Where A, of Eq (26) is expressed in Eq (27). 


Aa i-es [1-7] ol e) 


Where Q and Qs are the heat addition into the GT combustion 
chamber and the complementary firing process (the 
reheater) Figure 1. 
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Figure 8. Effect of primary zone temperature on 
environmental emissions 
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Figure 9. Effect of gasifier temperature and mass flow on the 
unit cost of electricity 


1000 


2500 


UCOE (N/kWh) at diff. gasifier mass flow 


23.0 240000 
—o—UCOE 
—s— Kalina Turbine PEC 220000 
22.5 
200000 € 
o 
zm 
A 
E 
22.0 180000 -£ 
5 
= 
160000 À 
z 
21.5 
140000 
21.0 120000 
0.25 0.3 0.35 0.4 


Ammonia mass fraction in water 
Figure 10. Effect of ammonia mass fraction on the unit cost of 
electricity 


1 


FI Abam et al. /Future Technology 


February 2024]| Volume 03 | Issue 01 | Pages 40-55 


Gasifier mass flow rate (kg/s) 


E 2 23 
5.350*10 


—t+— Total cost at varying gasifier temp. 
—s— Total cost at varying gasifier mass flow 


5.300*10° 


5.250x10* 


5.200x10* 


5.150x10* 


System total cost ($) (from biomass gas. temp. 


5.100x10* 


500 600 


3 


700 


3.5 4 4.5 5 " 
5.35510 


5.350x10* 


5.345x10° 


5.340x10* 


5.335*10° 


6 


5.330*10 


6 


5.325*10 


System total cost ($) (from biomass mass variation) 


5.320*10° 


800 900 1000 


Gasifier temperature (°C) 


Figure 11. Effect of gasifier temperature and mass flow on system total cost 


Table 13. Components of optimization functions and parameters 


Performance Optimization function Decision variables Optimization constraints 
index 
T. 
Wirc mo —-[(5). - 1] Tuc Tp 0.80 < LPC < 0.90 
PC 
T: 
Wupc rise ——[(r) - 1] Mpc Tp 0.80 < LPC < 0.90 
PC 
W, : 1 p 0.80 < LPC < 0.9 
eT Tha 0Cyl oN pr |1 — | Nec» Tp, 
(5) 8x7, < 16,1150 < Ty) < 1250 
W, ; 1 h 0.80 < LPC < 0.9,8 < f, < 
HPT 1it7CyT Mur |1 — = Neco ly, p 
Tp) 16, 1150 < T; < 1250 K 
Wkarrurb Mrg|h2g — hao] + Pog, pressure at Tg, hg 435 € Tj, < 439K 
[rios x M29]. [h29 Ez has] 
Qevp1 Tnzolhz1 — hz2l Tao Tsx -3 and -1.5 °C 
VPyenerator Tha 2| h45 = hys| Ti; 270< Teo < 282 
Qevp2 faga; — hao] Tao T44, -5.8 and -1.5 °C 
Cos, Total Cruet + Cony. + Y Le Described parameters 


k 


The total cost rate of the system and the cost rate due to 
environmental impact are presented in Eqs (28) and (29). The 
components of the optimization function, objective functions, 
constraints and key performance indices are depicted in 
Table 13. 


Crot = Cruel + Log F Yule (28) 


Cenv = Cno, ™no, + Cco,Mco, + CcoMco (29) 


Where the unit damage costs Cyo,, Cco, and Cco are taken as 
0.02186 $/kg, 6.863$/kg and 0.023 $/kg, respectively. 


The GA (genetic algorithm) was applied in the optimization 
due to its flexibility in simplifying multi-variable and multi- 
objective problems. In this analysis, Ninety groups of the 
Pareto-frontiers from the GA were described following the 
OBFs and the equivalent constraints. The 17th Pareto-front 
correspond to the optimum (y) and minimum cost rate 
calculated at 45.32 96 and 125.84 $/hr, separately. The 
corresponding parameters occur at a compression ratio (CR) 
of 8, with isentropic efficiencies of LPC, HPC, LPT, and HPT 
existing at 8896. Also, the intercooler exit temperature and the 
inlet temperature to the combustion chamber and reheater 
were calculated at 386.7K, 1140 K and 1240.3 K. The optimal 
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inlet pressures of the Kalina system were obtained at 20 kPa. 
The specific emissions rate at the optimum operating 
conditions are calculated at 123.34, 2.87E-07 and 0.214 kg/ 
MWh for CO2, CO and NOx, respectively. The improvement 
potential of the environmental discharges is for CO2, 1.56 96, 
NOx, 4.33 96 and CO, 3.65 96, with optimum hydrogen 
production. 


5. Conclusion 

A Kalina-based multigeneration cycle with increased 
heat addition to the PEM electrolyzer for hydrogen and water 
production was performed. The results are summarised as 
follows: The multigeneration energy system has a net energy 
and exergy efficiency of 53.48 and 50.05 96, respectively. 
Additionally, the system has improved topping cycle 
efficiency by 19.03 and 2.58 96, respectively, from the original 
values of 34.45 and 47.47 96. Atthe cycle level, the gas turbine 
contributed up to 85.81 96 of the ED. The combustion 
chamber and reheater contributed about 54.65 and 22.57 96 
of the total ED destructions. The system overall exergo- 
thermal index (ETI) stood at 1.713 with the gas turbine 
system alone the ETI was obtained as 2.106, which suggests a 
higher thermal impact on the environment without the 
addition of the lower cycles. The system has a hydrogen 
production rate of 0.1524 kg/h with cooling rate and COP 
calculated at 3498 kW and 4.304 respectively. The ESI for the 
system was relatively low, with a value of 0.5242. The low ESI 
arises from the large ED in the topping cycle, which 
contributed about 66 % of total system ED. The 
exergoeconomic result indicates a low average cost of energy 
from the gas turbine (0.836 $/GJ) compared to the Kalina 
subsystem (6.53 $/GJ) and vapour absorption system (3.24 
$/GJ). 
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Appendix 


Table 1. Summary of the energy and exergy balances, as well as the exergy of fuel and product for the system 
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Component Energy Balance Exergy Balance Exergy of Fuel Exergy of Product 
GT LPC mah, + Wipc = mih; Ey + Wipe = Ep + Epipc W,pc É, — Ey 

GT intercooler ^ rhjh; + Tü4gh,g = Mh3hz + Tiygl4o È, + Ès = È; + È, + Ep, ir È, — E, E45 — Egg 

GT HPC Thshs  Wypc = mh, Éz + Wiec = Ë, + Épupc Wurc É, — E, 

GT Hex Th, + Thy hy, = Tshs + mazhi2 É, + Ey, = Ès + Èz Enernex È — É È; — È, 

GT CC Me|LHV,| + mghs = Th;h; Es + Ès = É, + Eperee Es + Ès È, 

GT HPT titzh, = ghg + Wypr È, = Ég + Wypc + Wipe +Èp,upr E, — Es Wupc + Wipc 
GT REH Mg|LHVo| + ghg = Moho Ès + Ey = Eyo+ Épr ren Eg + Eo Eio 

GT LPT Tniohio = Mihi + Wier Ey = È + Wher + Éppr Èo — Ey Wher 

Kal. Vap Gen. Myghy, + M37hg7 = Myghy3 + Myzhy7 Ey, + B37 = Ey3 + E17 Épyapcen Ey, — Ei Ey7 — E37 
Kal. Sep 1 Mizhy7 = Mighyg + mazh22 Ey; = Exg + Ej; + Epxaisi Ey; Eis + Eo, 
Kal. Valve 1 Ma2hz2 = My3hz3 Ey, = Èn + Ep, Kal V1 E? E55 

Kal. Valve 2 mh; = mjhj, Ey, = Èn + Ep, Kal V2 Èz Ej, 

Kal. Turb. Mightg = zoho + MigWrga + M34Wp, + ra3 We, Éis = Éz + Éw, ea + Ew,, + Ewp, + Eben Éis — E29 Éw, i + Èw, + Ewp, 
Kal. Sep 2 Maohz9 = magh;g + mah; Eq = Éj + Ext + Epeatse Exo Ès + Eo, 
Kal. CND 1 Maghyg + Msohsg = Myghzq + Megheo Eng T És, = - Es + Foo + Épcupi Èz — Éjs Éso — Eso 
Kal. Valve 3 fh;gh;g = Maoh29 Ej, = = Es T Ep, Kal V3 Eng f É;s . 

Kal. Evap 1 Myghy, + mah. = Myghy3 mah Ey, + Ey = Ey3 + Ey7 + Epy.cen Ej; — E13 Ej; — E37 
Kal. HEX 1 My6h6 + Ms6hs6 = mhi; + Mz7hz7 Èz Ey = Ey + Ès7 + Èp kanxe Ézs — E, E57 — Ese 
Kal. HEX 2 M3zh32 + Mg5h35 = M33h33 + mucha Es * Ess = Ess nS E56 + Èb kai nex2 E32 — Eas E36 — Eas 
Kal. CND 2 M33h33 + Ms7hs7 = Mz4hzq + Msghsg B33 + Eg; = E34 + Esg + Ep rar cnp2 E33 — Esq Esg = Es; 
Kal. Pump 1 Tha 4ha4 + M34Wp, = Thashas E54 + Ew, = Es Epp, Ew, Eas — E34 
VAS desorber Mashys + Myghyg = magh;g + maghag + mashas Es + E44 = Ézs + Egg + Eas + Eppese Ézs — Éjg + Eng Ézs — Eas 
VAS HEX 3 Maghyz + Mashys = Maghgg + Myphae Ej zh Eis = Ej + Exe + Epvasuexs É,s — Eas Ej, — Eis 
VAS Valve 5 Maghyg = mazha; Bag = Ea; Epyasvs Es m ; 

VAS Pump 2 Tha5ha5 + M42Wp, = Mazhy3 Ej; Ew, = Ess Epp, Ewp, E43 — Ea 
VAS Absorber mua ha, omazha; + Meshes = Mazhy2 + Meshes Ea + Ej; + Éss = Ei; * Egg + Epasse Eq, +È Eee — Eos + Ez 
VAS EVP 2 Maghgo + Mezhez = Margi mg hg, Ezo + Egg = Es, + Eg. + Ep,vas ever Eo3 — [m Eq Ex 
VAS Valve 4 Mgghzq = Maoh4o E39 = Ego + Ep, VAS V4 E39 Eso 

VAS CND 3 Maghag + Meihe1 = Maghzq + Mezhe2 Eg + Eo, = E39 + Eo2 Épyascups Ézs — Egg Ész — Ee 
Water heater mash + Myshys = mih + Michie Èis + Èis = Ë + Èis + Epnearer Ey3 — Ex, Éis — Éis 
PEM Electr. Tszħs2 + Myghyg = Mszhsz3 + Ms4hsq Es2 + Big = Es3 + Esa + Ép giectrolyser Ész + E19 És + Esa 
PEM HEX Msghsg + Meoheo + Msohso = Msihs1 + Megheo Esg + Boo + Eso = Ës1 + Eoo + Ep,pem nex Esg + Éso — Eos És, — Eso 
PEM 02 Msahs4 = Msshss + Msehso E54 = Ess + Eso + Époxysep Esa Ess + Ese 
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